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ABSTRAGT 

Merging galaxy clusters leave long-lasting signatures on the baryonic and non-baryonic cluster con¬ 
stituents, including shock fronts, cold fronts. X-ray substructure, radio halos, and offsets between 
the dark matter and the gas components. Using observations from Chandra^ the Jansky Very Large 
Array, the Giant Metrewave Radio Telescope, and the Hubble Space Telescope, we present a mul¬ 
tiwavelength analysis of the merging Frontier Fields cluster MAGS J0416.1-2403 {z = 0.396), which 
consists of a NE and a SW subclusters whose cores are separated on the sky by ^ 250 kpc. We 
find that the NE subcluster has a compact core and hosts an X-ray cavity, yet it is not a cool core. 
Approximately 450 kpc south-south west of the SW subcluster, we detect a density discontinuity 
that corresponds to a compression factor of ^ 1.5. The discontinuity was most likely caused by the 
interaction of the SW subcluster with a less massive structure detected in the lensing maps SW of 
the subcluster’s center. For both the NE and the SW subclusters, the dark matter and the gas com¬ 
ponents are well-aligned, suggesting that MACS J0416.1-2403 is a pre-merging system. The cluster 
also hosts a radio halo, which is unusual for a pre-merging system. The halo has a 1.4 GHz power of 
(1.06 ±0.09) X 10^^ W Hz which is somewhat lower than expected based on the X-ray luminosity 
of the cluster. We suggest that we are either witnessing the birth of a radio halo, or have discovered 
a rare ultra-steep spectrum halo. 

Subject headings: G3lax.ies: clusters: individual: MACS J0416.1-2403 — Galaxies: clusters: intraclus¬ 
ter medium — X-rays: galaxies: clusters 


1. INTRODUCTION 

Galaxy clusters grow by merging with other clusters 
and by accreting smaller mass structures from the in- 
tergalactic medium. Signs of these interactions are im¬ 
printed in the ICM and detected in X-ray observations as 
cold fronts, shock fronts, turbulence, and ram pressure- 
stripped gas (e.g., Markevitch & Vikhlinin 2007; Randall 
et al. 2008a; Zhuravleva et al. 2015). Other footprints of 
cluster interactions can be seen in the radio band as halos 
and relics (e.g., Feretti et al. 2012). If the merger is not 
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in the plane of the sky, mergers can be detected in optical 
observations based on multiple peaks in the radial veloc¬ 
ity distribution and in the spatial galaxy distribution. 
Furthermore, comparison of X-ray and optical/lensing 
data also reveals signatures of merging events, most no¬ 
tably as offsets between the dark matter (DM) and the 
gas components of the colliding clusters (e.g., Marke¬ 
vitch et al. 2004; Glowe et al. 2006; Randall et al. 2008b; 
Merten et al. 2011; Dawson et al. 2012). 

Radio halos, which have been detected in some merg¬ 
ers, are diffuse synchrotron-emitting sources with low 
surface brightness, steep spectral indices (g < — 1, Jv oc 
z/*^), and typical sizes of ^ 1 Mpc. Because of their steep 
radio spectra, low-frequency radio observations play an 
important role in characterizing their properties. Two 
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main models have been proposed to explain the origin of 
such halos: 

• Reacceleration by turbulence: Large-scale turbu¬ 
lence generated during the merger event sup¬ 
plies the energy required to reaccelerate fossil cos¬ 
mic rays (CRs) back to relativistic energies, at 
which time they become synchrotron-bright (e.g., 
Brunetti et al. 2001; Petrosian 2001). 

• Acceleration by hadronic collisions: Inelastic colli¬ 
sions between CR protons trapped in the gravita¬ 
tional potential of a cluster (e.g. originating from 
supernova explosions, active galactic nuclei [AGN] 
outbursts, previous merger events, etc.) and ther¬ 
mal ICM protons give rise to a secondary popula¬ 
tion of CR electrons, which consequently are visi¬ 
ble at radio frequencies (e.g., Blasi & Colafrancesco 
1999; Dolag & EnBlin 2000). 

However, hybrid models have also been postulated, in 
which radio halos may be produced by turbulent reac¬ 
celeration of secondary particles resulting from hadronic 
proton-proton collisions (Brunetti & Blasi 2005; Brunetti 
& Lazarian 2011). The radio halos predicted by hybrid 
models are expected to be found in more relaxed clusters 
and to be underluminous for the masses of the hosting 
clusters (see also Brunetti & Jones 2014, for a review). 

Hadronic radio halo models predict that magnetic 
fields should either be different in clusters with and with¬ 
out halos, or that gamma-ray emission should be de¬ 
tected in clusters hosting halos (e.g., Jeltema & Profumo 
2011). However, there is no evidence that magnetic fields 
are stronger in clusters with radio halos (e.g., Bonafede 
2010), and there has been no conclusive gamma-ray 
detection with the Fermi telescope (Ackermann et al. 
2010). Therefore, current observational evidence disfa¬ 
vors hadronic models. 

A textbook example of a merging cluster with a bright 
radio halo is the famous Bullet cluster, IE 0657-56 (Elvis 
et al. 1992; Liang et al. 2000a; Markevitch et al. 2002; 
Shimwell et al. 2014). Chandra has revealed that this 
system has a bullet-like gas cloud moving through the 
disturbed ICM of the main cluster (Markevitch et al. 
2002). The gas cloud is the surviving core of a sub¬ 
cluster whose outer layers were ram pressure-stripped 
in the collision. Immediately ahead of the core there 
is a density discontinuity associated with a cold front, 
while further out in the same direction there is a den¬ 
sity discontinuity associated with a shock front (Marke¬ 
vitch et al. 2002; Owers et al. 2009). Recently, another 
shock front has been found on the opposite side of the 
Bullet cluster, behind the bullet-like gas cloud (Shimwell 
et al. 2015). Chandra observations of the Bullet clus¬ 
ter were followed up with optical/lensing data acquired, 
most notably, with the Very Large Telescope (VLT) and 
the Hubble Space Telescope (HST). Chung et al. (2010) 
have shown that the redshift distribution of the Bullet 
cluster is bi-modal, with two redshift peaks at z ^ 0.21 
and 2 : ^ 0.35, which imply a velocity difference ^ 3000 
km s“^ between the two merging subclusters. The gravi¬ 
tational leasing analysis has shown that the Bullet cluster 
is a dissociative merger, in which the essentially collision- 
less DM (ctdm < 0.7 g cm“^, Randall et al. 2008b) has 
decoupled from the collisional ICM (Markevitch et al. 


2004; Randall et al. 2008b). Lage & Farrar (2014) have 
combined the X-ray and leasing data with numerical sim¬ 
ulations to constrain the merger scenario of the Bullet 
cluster. 

Here, we present results from Chandra^ Jansky Very 
Large Array (JVLA), and Giant Metrewave Radio Tele¬ 
scope (GMRT) observations of another merging galaxy 
cluster: the HST Frontier Fields^ cluster MACS J0416.1- 
2403 (z=0.396, Ebeling et al. 2001; Mann & Ebeling 
2012). Optical and leasing studies of this cluster have 
been presented recently by, e.g., Zitrin et al. (2013); 
Schirmer et al. (2014); Jauzac et al. (2014); Zitrin et al. 
(2015); Grillo et al. (2015). Most recently, Jauzac et al. 
(2015) mapped the DM distribution of this merging sys¬ 
tem using strong- and weak-lensing data. Their analysis 
revealed two mass concentrations associated with the two 
main sub clusters involved in the merger, plus two smaller 
X-ray-dark mass structures NE and SW of the cluster 
center. Jauzac et al. (2015) combined their lensing data 
with archival Chandra observations to study the offsets 
between the DM and the gas components. They report 
good DM-gas alignment for the NE subcluster, but a sig¬ 
nificant offset for the SW one. Based on the lensing and 
X-ray results, Jauzac et al. (2015) proposed two possible 
scenarios for the merger event in MACS J0416.1-2403 - 
one pre-merging and one post-merging - but were unable 
to distinguish between the two. 

The analysis presented herein combines significantly 
deeper Chandra observations with recently acquired 
JVLA and GMRT data and with the optical/lensing re¬ 
sults reported by Jauzac et al. (2015) to improve our un¬ 
derstanding of the merger event in MACS J0416.1-2403. 
The Chandra JVLA, and GMRT observations and data 
reduction procedures are described in Section 2. In Sec¬ 
tion 3 we discuss the background analysis of the Chandra 
data, and in Sections 4, 5, 6, and 7 we present the X-ray 
results. The radio results are presented in Section 8. In 
Section 9 we use the deeper Chandra data to revise the 
DM-gas offsets reported by Jauzac et al. (2015). In Sec¬ 
tion 10 we discuss the implications of our findings for the 
merger scenario of MACS J0416.1-2403. In Section 11 we 
examine the origin of the radio halo. A summary of our 
results is provided in Section 12. 

Throughout the paper we assume a ACDM cosmology 
with Hq = 70 km s~^ Mpc“^, Dm = 0.3, and Da = 0.7. 
At the redshift of the cluster, 1' corresponds to approx¬ 
imately 320 kpc. Unless stated otherwise, the errors are 
quoted at the 90% confidence level. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Chandra 

The HST Frontier Fields cluster MACS J0416.1-2403 
was observed with Chandra for 324 ks between June 2009 
and December 2014. A summary of the observations is 
given in Table 1. 

The data were reduced using Cl AO v4.7 with the cali¬ 
bration files in CALDB v4.6.5. The particle background 
level of the VFAINT observations was lowered by filter¬ 
ing out cosmic ray events associated with significant fiux 
in the 16 border pixels of the 5x5 event islands. Soft 
proton flares were screened out from all observations us¬ 
ing the deflare routine, which is based on the Ic.clean 

® http://www.stsci.edu/hst/campaigns/frontier-fields/ 
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code written by M. Markevitch. The flare screening was 
done using an off-cluster region at the edges of the field 
of view (FOV), which was selected to include only pixels 
located at distances greater than 2.5 Mpc from the clus¬ 
ter center and to exclude any point sources detectable 
by eye. The clean exposure times following this cleaning 
are listed in Table 1. 

The inspection of an ObsID 10446 spectrum extracted 
from an off-cluster region showed a significant high- 
energy tail, which indicates that this observation is con¬ 
taminated by flares that were not detected by the fil¬ 
tering routine. Given the relatively short exposure time 
of this observation, we decided to exclude ObsID 10446 
from our analysis rather than attempt to model the flare 
component. 

We reprojected the other five observations to a com¬ 
mon reference frame, and created merged images in the 
energy bands 0.5 — 2, 0.5 — 3, 0.5 — 4, 0.5 — 7 keV, and 
2 — 7 keV. The exposure-corrected, vignetting-corrected 
0.5 — 4 keV mosaic image is shown in Figure 1. To de¬ 
tect point sources, we also created merged exposure map- 
weighted PSF images (EOF = 90%) in each of the five 
energy bands. Point sources were detected in the indi¬ 
vidual bands with the CIAO task wavdetect, using the 
merged maps, wavelet scales ofl, 2, 4, 8, 16, and 32 pix¬ 
els, a sigma threshold of 5 x 10“^, and ellipses with 5cr 
axes. A few additional point sources that were missed 
by wavdetect were selected by eye and excluded from 
the data. All point sources detected by wavdetect were 
excluded from the data very conservatively, using the el¬ 
liptical region that covered the largest area among the 
five elliptical regions identified for the different energy 
bands. 

2.2. JVLA 

JVLA observations of MACS J0416.1-2403 were ob¬ 
tained in the L-band in the BnA, CnB, and DnC ar¬ 
ray configurations. The observations were recorded with 
the default wide-band L-band setup, giving 16 spectral 
windows each having 64 channels, covering the entire 1- 
2 GHz band. An overview of the observations is given in 
Table 2. 

The data were calibrated with the Common Astron¬ 
omy Software Applications^ (CASA) package version 
4.2.1. As a first step we flagged data affected by radio 
frequency interference (RFI) using AOFlagger (Offringa 
et al. 2010), after correcting for the bandpass shape. 
Data affected by other problems and antenna shadowing 
were flagged as well. After flagging, we applied the pre¬ 
determined elevation dependent gain tables and antenna 
offset positions. The data were Hanning smoothed. 

We determined initial gain solutions using 10 chan¬ 
nels at the center of the spectral windows for the 
primary calibrators 3C147 and 3C138. We then re¬ 
calibrated the bandpass and obtained delay terms 
(gaintype=^K^ ) using the unpolarized calibrator 3C147. 
A next step consisted of calibration of the cross-band de¬ 
lays (gaintype= ^KCROSS ^) using the polarized calibrator 
3C138. We then solved again for the gains on the primary 
calibrators but now using all channels. The channel de¬ 
pendent leakage corrections were found using 3C147 and 
the polarization angles were set using 3C138. The gains 

^ http://casa.nrao.edu 


were again re-determined for all calibrator sources which 
included the phase calibrator J0416-1851. Finally, the 
flux density-scale was bootstrapped from the primary 
calibrators to J0416-1851 and the calibration solutions 
were applied to the target field. 

To refine the calibration for the target field, we per¬ 
formed three rounds of phase-only self-calibration and 
two final rounds of amplitude and phase self-calibration. 
W-projection (Cornwell et al. 2008, 2005) was employed 
during the imaging, taking the non-coplanar nature of 
the array into account. The self-calibration was indepen¬ 
dently performed for the three different datasets from the 
different array configurations. The full bandwidth was 
imaged using MS-MFS clean (nterms=2; Ran & Corn- 
well 2011). Briggs (1995) weighting with a robust fac¬ 
tor of -0.75 was used for self-calibration. Clean masks 
were used, which were made with the PyBDSM^ source 
detection package. 

After the self-calibration, the three datasets were com¬ 
bined. One final round of self-calibration was carried 
out on the combined dataset. The final images were cor¬ 
rected for the primary beam attenuation. 

In addition, we imaged the dataset using an inner uv- 
range cut of 4.3 kA (corresponding to a scale of about 1', 
or 320 kpc). This model was then subtracted from the 
visibility data to allow a search for diffuse emission in the 
cluster. We imaged the dataset with the emission from 
compact sources subtracted using a Gaussian uv-taper 
of 20" and employing multi-scale clean (Cornwell 2008). 
An overview of the resulting image properties, such as 
root-mean-square (rms) noise and resolution are given in 
Table 3. 

2.3. GMRT 

GMRT 610 MHz data for MACS J0416.1-2403 were 
collected on December 3, 2013 using 26 antennas. The 
on-source time was 5.0 hrs, and a total bandwidth of 
33.3 MHz (RR correlations only), split into 512 chan¬ 
nels, was recorded. NRAO’s Astronomical Image Pro¬ 
cessing System (AIPS) was used to carry out the initial 
calibration of the visibility dataset. The primary cali¬ 
brator 3C 147 was used to set the flux density scale and 
derive the bandpass solutions for all the antennas. The 
source J0409-179 was observed for 6-mins scans at inter¬ 
vals of ^ 25 mins and used as the secondary phase and 
gain calibrator. About 20% of the data, mostly from 
short baselines, were affected by RFI and subsequently 
flagged. Gain solutions were obtained for the calibrator 
sources and together with the bandpass solutions applied 
to the target field. In total, 480 of the 512 channels were 
used; the rest of the channels were discarded as they were 
too noisy due to the bandpass roll-off. The 480 channels 
were averaged down to 48 channels to reduce the size of 
the data. 

The visibility data were then imported into CASA to 
refine the calibration via the process of self-calibration. 
Three rounds of phase-only self-calibration and two final 
rounds of amplitude and phase self-calibration were ap¬ 
plied. The phase-only self-calibration was carried out 
on a 30 s timescale. The amplitude and phase self¬ 
calibration was carried out on a 2 min timescale, pre- 

® http://dl.dropboxusercontent.com/u/1948170/html/index. 
html 
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Table 1 

Chandra observations 


ObsID 

Observing mode 

CCDs on 

Starting date 

Tot^ yme 

Cle^g jime 

10446 

VFAINT 

0,1, 2, 3, 6 

06-07-2009 

15.8 

15.8 

16236 

VFAINT 

0,1, 2, 3, 6 

08-31-2014 

39.9 

38.6 

16237 

FAINT 

0,1, 2, 3, 6 

11-20-2013 

36.6 

36.3 

16304 

VFAINT 

0,1, 2, 3, 6 

06-10-2014 

97.8 

95.2 

17313 

VFAINT 

0,1, 2, 3, 6 

11-28-2014 

62.8 

59.9 

16523 

VFAINT 

0,1, 2, 3, 6 

12-19-2014 

71.1 

69.1 


Table 2 

JVLA L-band Observations 



Bn A-array 

CnB-array 

DnC-array 

Observation dates 

25 Jan, 2014 

12 Jan, 2015 

Sep 19, 2014 

Frequency coverage (GHz) 

1-2 

1-2 

1-2 

On source time (hr) 

0.6 

3.2 

1.5 

Correlations 

full stokes 

full stokes 

full stokes 

Channel width (MHz) 

1 

1 

1 

Visibility integration time (s) 

2 

3 

5 


Table 3 

Radio imaging parameters 


Instrument 

Weighting 

Resolution 
(arcsec x arcsec) 

r.m.s. noise 
(/xJy beam“^) 

JVLA 1.25 GHz 

Briggs -0.75 

7.8 X 5.5 

14 

JVLA 1.25 GHz 

uniform + 20" Gaussian taper 

18 X 18 

20 

GMRT 610 MHz 

Briggs 0.5 

7.6 X 4.0 

54 

GMRT 610 MHz 

uniform + 20" Gaussian taper 

20 X 20 

0.36 X 10^ 

VLITE 340 MHz 

Briggs 0.0 

45.7 X 31.1 

2.1 X 10^ 


applying the phase-only solutions before solving for the 
amplitude and phases on the longer 2 min timescale. 
We employed W-projection during the imaging and used 
nterms=2. 

A map with Briggs (robust=0.5) weighting of the field 
was produced, which resulted in a resolution of 7.6" x 4.0" 
and an rms noise level of 54 /iJy beam“^ (Table 3). Simi¬ 
larly to the JVLA L-band data reduction, we imaged the 
dataset using an inner uv-range cut of 4.3 kA to obtain 
a model of the compact sources in the field. After sub¬ 
tracing this model from the uv-data, we re-imaged with 
uniform weighting and a Gaussian uv-taper to obtain a 
beam size of 20". 

2.4. VLITE 

A new commensal observing system called the 
VLA Low Band Ionospheric and Transient Experiment 
(VLITE)^ has been developed for the NRAO JVLA 
(Clarke et ah, in prep.) and was operational during the 
observations of MACS J0416.1-2403. The VLITE cor¬ 
relator is a custom designed DiEX software correlator 
(Deller et al. 2011). The system processes 64 MHz of 
bandwidth centered on 352 MHz with two second tem¬ 
poral resolution and 100 kHz spectral resolution. VLITE 
operates during nearly all pointed JVLA observations 
with primary science goals at frequencies above 1 GHz, 

® http://vlite.nrao.edu/ 


providing data simultaneously for 10 JVLA antennas us¬ 
ing the low band receiver system (Clarke et al. 2011). 

We processed the CnB configuration VLITE observa¬ 
tion of MACS J0416.1-2403 from 12 January 2015 using a 
combination of the Obit software package (Cotton 2008) 
and AIPS (van Moorsel et al. 1996). VLITE data at fre¬ 
quencies u >360 MHz were removed due to the presence 
of strong REI from a satellite downlink that is present 
during most operations. The data at u <360 MHz were 
flagged using the AIPS program RELAG to remove the 
majority of the remaining REI. The bandpass was flat¬ 
tened using observations of several calibrators taken near 
the time of the observations (3C48, 3C138, 3C147, and 
3C286). These calibrator sources were taken in the same 
primary observing band as MACS J0416.1-2403. The de¬ 
lays were determined using these same calibrators. The 
delay corrected data then were flux calibrated using these 
same calibrator sources. No phase calibrator is required 
for low band observations due to the large field of view 
(EWHM ^ 2.3° at 320 MHz) and typical presence of 
sufficiently strong sources in the field that allow for self¬ 
calibration. 

The initial imaging of the target field revealed addi¬ 
tional REI that was identified in the residual data set 
after all compact sources had been subtracted from the 
uv data. These baselines were excised from the full data 
set and we further refined the calibration through four 
rounds of phase-only self-calibration. Non-coplanar ef- 
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Figure 1. Logarithmic-scaled Chandra surface brightness map of the HST Frontier Fields Cluster MACS J0416.1-2403, in the energy 
band 0.5 — 4 keV. The image was binned by a factor of 4 (1 pixel 2"), exposure-corrected, vignetting-corrected, and smoothed with a 
2D Gaussian kernel of width 1 pixel xl pixel. The large image shows the mosaic of all the ObsIDs summarized in Table 1, while the top 
image zooms in on the brightest X-ray regions of the cluster, in a region of size 2.9 X 2.9 Mpc^. The dashed grey circle marks the boundary 
of the -R 500 region of the cluster. 







6 


Ogrean et al. 


fects were taken into account in the imaging steps us¬ 
ing small facets to make a fly-eye image of the full 
field out to the FWHM with additional facets placed on 
bright sources out to 20° from the field center. Data 
were imaged using small clean masks placed around all 
sources. The final VLITE image of MACS J0416.1-2403 
has a resolution of 45.7" x 31.1" and an rms noise of 
2.1 mJy beam“^. The image was made using a Briggs 
robust parameter of 0 (see Table 3). 

3. X-RAY FOREGROUND & BACKGROUND MODELLING 

We analyzed the Chandra data using the Group F 
stowed background event files, which are appropriate for 
observations taken after September 21, 2009. The parti¬ 
cle background was subtracted from the data, while the 
foreground and background sky components were mod¬ 
eled from off-cluster regions. For each of the cluster 
ObsIDs, we created corresponding stowed background 
event files, applied the VFAINT cleaning to those as¬ 
sociated with cluster observations taken in this observ¬ 
ing mode, and renormalized the stowed background ob¬ 
servations such that the count rates in the energy band 
10 — 12 keV matched those of the corresponding obser¬ 
vation in the same energy band. We modeled the sky 
components as the sum of unabsorbed thermal emission 
from the Local Hot Bubble (LHB; APEC component 
with T ^ 0.1 keV; Snowden 1998), absorbed thermal 
emission from the Galactic Halo (GH; APEC component 
with T ^ 0.2 keV; Henley et al. 2010), and absorbed 
non-thermal emission from undetected point sources in 
the FOV (power-law component with index 1.41 ± 0.06; 
De Luca & Molendi 2004). The redshifts of the fore¬ 
ground components were set to 0, while the abundances 
were set to solar values assuming the abundance table of 
Anders & Grevesse (1989). The hydrogen column den¬ 
sity was fixed to 2.89 x 10^^ atoms cm“^ (Kalberla et al. 
2005). The Chandra spectra were fit in the energy band 
0.5 — 7 keV using XSPEC vl2.8.2 - the lower limit was 
chosen to avoid calibration uncertainties at low energies, 
and the upper limit to increase the signal-to-noise ratio 
(SNR). 

When fitting the Chandra spectra in the energy band 
0.5 — 7 keV, the ^0.1 keV LHB component cannot be 
constrained. Instead, we constrain this component using 
a ROentgen SATellite (ROSAT) All-Sky Survey (RASS) 
spectrum^^ corresponding to an annulus with radii 0.15 
and 1.0 degrees around the cluster. The ROSAT spec¬ 
trum was fit in the energy band 0.1 — 2.4 keV, with the 
normalization of the power-law background component 
fixed to 8.85 x 10“^ photons keV“^ cm“^ s“^ arcrnin”^ 
(Moretti et al. 2003). The temperatures and normaliza¬ 
tions of the LHB and GH components were free in the 
fit. The results of the fit to the ROSAT spectrum are 
presented in Table 4. 

Based on the ROSAT results, the parameters of the 
LHB components were fixed for the Chandra spectra. 
The five Chandra spectra were then fitted in parallel, as¬ 
suming that they are described by the same foreground 
model. The power-law normalizations, on the other 
hand, were left to vary independently, in order to ac¬ 
count for the varying exposure time across the merged 

https://heasarc.gsfc.nasa.gov/cgi-bin/Tools/xraybg/ 
xraybg.pi 


observation that was used for point source detection. 
The energy sub-band 0.5 — 0.8 keV was ignored for Ob- 
sID 16237 due to negative spectral residuals caused by a 
change in the spectral shape of the background compo¬ 
nent that is not filtered by the VFAINT cleaning, which 
occurred between 2009 (the date of the Group F back¬ 
ground files) and 2013 (the date of the observation; A. 
Vikhlinin, priv. comm.). The subtraction of the instru¬ 
mental background resulted in small line residuals near 
the spectral positions of the Al Ka {E ^ 1.49 keV), Si Ka 
{E « 1.75 keV), and An M(a,/3 {E ^ 2.1 keV) fluores¬ 
cent instrumental lines. Therefore, we excluded from the 
spectra very narrow bands {AE = 0.10 keV) surrounding 
these lines. 

The best-fitting foreground and background param¬ 
eters are summarized in Table 4. The spectra were 
grouped to have at least 1 count per bin, and the fit 
was done using the extended G-statistic^^ (Gash 1979). 

4. GLOBAL X-RAY PROPERTIES 

MAGS J0416.1-2403 has been previously identified as 
a merging galaxy cluster at 2 : = 0.396 (Mann & Ebeling 
2012). The brightness of the NE subcluster core is sig¬ 
nificantly more peaked than that of the SW subcluster, 
which instead has a rather flat brightness distribution 
(Figure 1). 

We determined the average cluster temperature by ex¬ 
tracting spectra from a circular region of radius 1.2 Mpc 
(approximately R5007 Sayers et al. 2013) centered at 
a = 4^16^08.8" and J = -24° 04'14.0" (J2000). The 
spectra extracted from the five Chandra ObsIDs were 
fit in parallel, assuming the cluster emission is described 
by a single-temperature thermal component^^. For this 
global fit, the sky background parameters were fixed 
to the values listed in Table 4, and the redshift was 
fixed to 0.396. We found T = lO.Odlg ^g keV and 

Lx, 0.1-2.4 keV = (9.14T0.10) X 10^"^ erg s“^. The average 
cluster parameters are shown in Table 5. 

5. MAPPING THE GLUSTER’S X-RAY PROPERTIES 
5.1. Mapping of the ICM Temperature 

We mapped the cluster properties using CONTBIN 
(Sanders 2006).^^ The merged Chandra image was adap¬ 
tively binned in regions of ^ 3600 source plus sky back¬ 
ground counts in the energy band 0.5 — 7 keV. The bins 
follow the surface brightness contours of the cluster, thus 
minimizing possible gas mixing in bins located near den¬ 
sity discontinuities. Instrumental background and to¬ 
tal spectra were extracted from regions corresponding to 
each bin, and grouped to have at least one count per 
bin. The net spectra were modelled as the sum of sky 
background components plus a single temperature ab¬ 
sorbed thermal model (phabsxAPEC) with free temper¬ 
ature and normalization. The metallicity was fixed to 
Z = 0.24Zq (see Section 4). The sky background emis¬ 
sion was fixed to the model summarized in Table 4. The 

^ ^ https://heasarc.gsf c.nasa.gov/xanadu/xspe c/manual/ 

XSappendixStatistics.html 

We tried adding an additional thermal component to fit the 
cluster emission, but the parameters of the second component could 
not be constrained. 

We also created temperature maps using the codes described 
by Ghurazov et al. (2003) and Randall et al. (2008a), and obtained 
consistent results at the 90% confidence level. 
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Table 4 

Foreground and background spectral models 




Chandra 




ROSAT 

Component 






LHB 

0.10 ±0.01 

C Q q9 + U.t )2 
».OZ_o 83 

X 

10 -'^^ 

0.10 ±0.01 

9.32+^;“ X 10-7 

GH 

0 

*^•^^-0.08 

9 qi=;+1-20 
^•'^^-1.66 

X 

10-7 

n OO+0.05 
0'-^^-0.04 

X 10-^ 

power-law 


4 20“^®'®® 

73 

X 

10-7 

- 

- 

ObsID 16236 




power-law 


9 1 7+0.76 
' -0.86 

X 

10-7 

- 

- 

ObsID 16237 




power-law 


9 4c;+0-60 

X 

10-7 

- 

- 

ObsID 16304 




power-law 


o 

1 + 

o o 

bo bo 

X 

10-7 

- 

- 

ObsID 17313 




power-law 


4 OO+0.71 
^.oo_o 69 

X 

10-7 

- 

- 

ObsID 16523 





^ Temperature, in units of keV. 

^ Normalization, in units of cm“^ arcmin“^ for the thermal components, and in 
units of photons keV“^ cm“^ s“^ arcmin“^ for the non-thermal components. 
Fixed parameter. 


Table 5 

Cluster properties in R^oo 

T (keV) 

z (z©) 0.24l°;“ 

i:^x.o.i-2.4keV (erg s-i) (9.14 ±0.10) x 10^-^ 

spectra extracted from the five Chandra datasets were 
fitted in parallel, with the temperatures and normaliza¬ 
tions linked between the different ObsIDs. The fits were 
done using the extended C-statistic. The resulting tem¬ 
perature map is shown in Figure 2. 

The temperature is high throughout the cluster, which 
causes the uncertainties on the measurements to be high. 
The temperature map does not show significant struc¬ 
ture. Within the uncertainties, the temperatures out to 
^ 400 kpc from the cluster center are consistent with 
- 10 keV. 

5.2. Mapping of the ICM Pressure and Entropy 

The electron number density can be derived either from 
the model normalization of the fitted spectrum, or from 
the surface brightness in the regions of interest. Extract¬ 
ing the number density from either of these quantities 
requires assumptions about the cluster geometry. Here, 
we estimated the electron number density from the sur¬ 
face brightness; we note that our results are unchanged 
when deriving the density from the spectral normaliza¬ 
tion. 

The surface brightness is essentially proportional to the 
emission measure, 

c cx y nidi (1) 

where the integration is done along the line of sight and 
He is the electron number density. More accurately, the 
surface brightness is also temperature-dependent. This 
dependence introduces an uncertainty of < 10% in the 
pressure and entropy maps, which is less than the un¬ 
certainties on the temperatures and does not affect our 
results. 

We approximated the density as and defined the 


pseudo-pressure and pseudo-entropy as: 

^’ = 2^Co%2keV> (2) 

K = TCoP2^,y. (3) 

The pseudo-pressure and pseudo-entropy maps are in¬ 
cluded in Figure 2. 

Strong jumps in pressure and entropy would indicate 
the presence of strong density discontinuities in the ICM. 
However, neither the pressure map, nor the entropy map 
of MACS J0416.1-2403 shows evidence of such strong 
jumps between the regions. We note that while a pres¬ 
sure and entropy discontinuity can be seen between re¬ 
gions #3 and #0, the size of region #0 is too large for 
this result to truly indicate that there is a density dis¬ 
continuity at the boundary of the two regions. 

In Table 6 we list the best-fitting temperatures, spec¬ 
tral normalizations, and 0.5 — 2 keV count rates for 
the regions in Figure 2. The bin numbers necessary 
to relate Table 6 and Figure 2 are shown in Fig¬ 
ure 3. An interactive version of the temperature map, 
which combines Figures 2, 3 and Table 6 is avail¬ 
able at http://hea-www.cfa.harvard.edu/~gogrean/ 
interactive/MACSJ0416_Tmap.html. 


6 . PROPERTIES OF THE INDIVIDUAL SUBCLUSTERS 

The merging history of the individual subclusters is 
closely related to their degree of disturbance. In the fol¬ 
lowing two sections, we perform the imaging analysis of 
the NE and SW subclusters in order to characterize their 
merging states. Eor both subclusters, we created surface 
brightness profiles in sectors centered on the respective 
X-ray peaks, and attempted to model the profiles with 
isothermal ^d-models (Cavaliere & Eusco-Eemiano 1976, 
1978): 


^(r) = 



( 4 ) 


where So is the central surface brightness, Tc is the core 
radius, and r is the radius from the cluster centre. 
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Temperature (keV) 

6 8 10 12 14 16 18 



14.40s 09.60s 4hl6m04.80s 

RA (J2000) 


Pressure (arbitrary units) 
0.0001 0.0025 0.0050 0.0075 



-24°03'00.0" 


06'00.0" 


14.40s 


09.60s 4hl6m04.80s 
RA (J2000) 


Entropy (arbitrary units) 
10 ^ 10 ^ 



14.40s 09.60s 4hl6m04.80s 

RA (J2000) 


Figure 2. Temperature, pseudo-pressure, aud pseudo-eutropy maps of MACS J0416.1-2403. lu all maps, black elliptical regious mark the 
regious from which poiut sources have beeu removed. 



Figure 3. Left: The map shows the associatiou of regious iu Figure 2 with the fit values iu Table 6. Numbers are regiou uumbers. Black 
ellipses mark the regious where poiut sources have beeu removed. Right: Same surface brightuess map as iu Figure 1, with overlaid regious 
used for mappiug the physical properties of the ICM. 


Table 6 

Best-fittiug temperatures, uormalizatious, aud 0.5 — 2 keV couut rates for the regious iu Figure 3. 


Biu uumber 


ATb 



5x" 


0 

' •»o_o 92 

(2.24 ±0.06) 

X 

lO-o 

(7.39 ± 0.07) X 

10 “® 

1 

10.7412-10 

(1.70 ±0.05) 

X 

10-1 

(3.38 ±0.06) X 

10 “® 

2 

13.2511“^ 

(2.31 ±0.07) 

X 

10-1 

(4.34 ± 0.08) X 

10 “® 

3 

10.32l(;9i 

(2.77 ±0.08) 

X 

10-1 

(5.30 ± 0.09) X 

10 “® 

4 

11.6912-10 

(1.63 ±0.05) 

X 

10-3 

(2.96 ± 0.06) X 

10-1 

5 

9.73llg 

(9.80 ±0.29) 

X 

10-1 

(1.84 ± 0.04) X 

10-1 

6 

15.3913;^ 

(9.03 ± 0.23) 

X 

10-1 

(1.62 ±0.03) X 

10-1 

7 


(7.13 ±0.22) 

X 

10-1 

(1.27 ±0.03) X 

10-1 

8 

12.18l2;« 

7.101(1;22 X 

10 “^ 

(1.31 ±0.03) X 

10-1 

9 

11 43 + 2-29 

(5.51 ±0.17) 

X 10-1 

(1.02 ±0.02) X 

10-1 

10 

8 41+1-37 

(4.12 ±0.13) 

X 10-1 

(8.02 ±0.15) X 

10 -® 


^ Temperature, iu uuits of keV. 

^ Spectral uormalizatiou, iu uuits of cm“^ arcmiu“^. 

0.5 — 2 keV couut rate, iu uuits of photous cm“^ s“^ arcmiu“^. 
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Simple /d-models provide good descriptions of the X- 
ray profiles of galaxy clusters that do not have strong 
ICM temperature gradients (e.g., Ettori 2000b,a). Based 
on Figure 2, there are indeed no strong tempera¬ 
ture gradients in the NE and the SW subclusters of 
MACS J0416.1-2403. 

Before fitting, all the profiles were binned to a mini¬ 
mum of 1 count/bin. The fits used Cash statistics. Bit¬ 
ting was done using a modified version of the PROEEIT 
package^^ (Eckert et al. 2011). All the surface bright¬ 
ness profiles presented in the following subsections are 
in the energy band 0.5 — 4 keV. The profiles were in¬ 
strumental background-subtracted, and exposure- and 
vignetting-corrected. 


6.1. NE Subduster 


The sector used to create the surface brightness profile 
of the NE subcluster and the fits to this profile are shown 
in Eigure 4. The sky background surface brightness level 
was determined by fitting a constant to the outer bins 
(radii 5.0' —10.7') of the profile. The sky background level 
was kept fixed in the following fits. The very central part 
of the profile is significantly underestimated by the best¬ 
fitting /3-model. Instead, the profile is described better 
by a double /3-model: 


Sir) = ^ 5o,i 


i=l,2 


1 + 1 — 

^c, i 


-3/3+0.5 


(5) 


Double /3-models provide good representations of 
cooling-core clusters (e.g., Jones & Eorman 1984; Ota 
et al. 2013). However, the core of the NE subcluster 
in MACS J0416.1-2403 is far from cool, having a tem¬ 
perature > 10 keV based on the temperature map in 
Eigure 2. To constrain a possible cooler component, we 
examined the probability that the gas in bin #4 (see 
Eigure 3) has two phases: one corresponding to a “hid¬ 
den” cool core, and another corresponding to the hot 
plasma that increases the average core temperature to 
>10 keV. We modelled the spectrum of bin 7^4 with a 
two temperature APEC model. The normalizations of 
the two components were free in the fit, as was one of 
the temperatures; the temperature of the second, cooler 
thermal component was fixed to 5 keV. The fit sets an 
upper limit of 2.82 x 10“^ cm“^ arcrnin”^ on the nor¬ 
malization of the cooler component (corresponding to a 
luminosity over 5 times lower than that of the NE sub¬ 
cluster), and provides no improvement in the statistics of 
the fit. A cool component with a temperature < 5 keV 
would need to be even fainter. Cooler gas in the NE 
core could be masked by inverse Compton (IC) emission 
from the active galactic nucleus (AGN) hosted by the 
NE subcluster brightest cluster galaxy (BCG). To exam¬ 
ine this possibility, we measured the temperature in an 
annulus with radii 32 and 64 kpc around the AGN in 
the NE BGG. The best-fitting temperature in this annu¬ 
lus is very high, 15.65l3;22 keV. In a smaller circle with 
a radius of 35 kpc around the NE core, the best-fitting 
temperature is IO. 54 I 211 keV, and consistent with the 
temperature calculated in an annulus around the core. 
Therefore, we find no evidence that the NE core is cool. 


Available upon request. 


Erom the best-fitting double ;d-model, we calculated 
the central density of the NE core to be (1.4 ± 0.3) x 
10“^ cm“^. The density and temperature of the NE core 
hence imply a cooling time of 3. 5 ^ 0.9 Gyr, which further 
supports our conclusion that the NE subcluster cannot 
be classified as a cool core cluster based on currently 
available X-ray data. 

To investigate if the double ;d-model shape of the NE 
profile is caused by substructure in a particular direc¬ 
tion, we divided the NE sector shown in Eigure 4 into 
three subsectors, and modeled the profile of each subsec¬ 
tor with [3- and double /d-models. The fits are shown in 
Eigure 5. Eor each of the subsectors, a double /3-model 
describes the profile better than a single /3-model at a 
confidence level > 99.99%. 


6.2. SW Subcluster 


The surface brightness profile of the SW subcluster and 
the sector used in the extraction of this profile are shown 
in Eigure 6. As for the NE profiles, the sky background 
profile was modelled by fitting a constant to the outer 
bins of the profile, in the range r = 4.5' — 10.7'. Unlike 
the NE profiles, the SW profile is modelled well by a 
simple ;d-model. The fit is shown in Eigure 6. 

A weak edge in the profile can be seen at r ^ 1.5'. 
We divided the profile into three subsectors with equal 
opening angles (60°) to check whether the edge is seen 
along all three directions, and found that it is present 
only in the central subsector (position angles 280° — 340°, 
measured from the W in a counterclockwise direction). 
To describe it, we fitted a projected broken power-law 
elliptical density model to the surface brightness profile 
between r = 0.5' — 4.0'. The density model is defined 
as: 


n{r) = < 




if r < rd 
if r > Td 


(6) 


where n is the electron number density, C is the density 
compression, and rd is the radius of the density jump. 
The fit is shown in Eigure 7. If the density discontinu¬ 
ity is a shock front, then its magnitude corresponds to 
a shock with Mach number Ai = 1.40lo!i2* Unfortu¬ 
nately, the count statistics are too poor to allow us to 
distinguish between a cold front and a shock front based 
on the temperature jump, and thus we cannot determine 
the nature of the surface brightness edge. 


7. SUBSTRUCTURE IN THE ICM 


The search for substructure in the IGM is motivated by 
the identification in the lensing maps of two less massive 
structures in addition to the main NE and SW subclus¬ 
ters (Jauzac et al. 2015). The positions of these mass 
structures are shown in Eigure 8 and denoted by SI and 
S2 for consistency with the notation of Jauzac et al. 
(2015). In the analysis of Jauzac et al. (2015), SI and S2 
were found to be X-ray-dark; however, the X-ray data 
presented here are ^ 6 times deeper, which would make 
it easier to observe ICM substructure. 

We searched for substructure using the unsharp- 
masked image of the cluster, which was created by divid¬ 
ing the difference of two 0.5 — 4 keV fluxed images con¬ 
volved with Gaussians of widths 4" and 10" by their sum. 
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R [arcmin] 



7x10“^ 10“^ 2x10“^ 1 2 3 4 5 

R [arcmin] 


Figure 4. Top: Sector used to model the surface brightness of 
the NE subcluster. Annuli are drawn only to guide the eye and 
do not reflect the actual bin size used for the surface brightness 
profiles. Middle: ^-model fit to the surface brightness profile of 
the NE subcluster. For clarity, the profile shown in this plot was 
binned to a uniform signal-to-noise ratio of 5. Bottom: Double 
^-model fit to the surface brightness profile of the NE subcluster. 
The individual /3-models are shown with dashed lines. For clarity, 
the profile shown in this plot was binned to a have 200 counts per 
bin. The best-fitting model parameters are listed in the middle and 
bottom plots; radii units are arcmin, and surface brightness units 
are photon cm“^ s“^ arcmin“^. Fixed parameters are shown with 
a superscripted f- The bottom panel shows the residuals of the fit. 


The resulting image, shown in Figure 8, highlights sub¬ 
structure on scales of ^ 20 — 50 kpc.^^ There is no excess 
X-ray emission at the positions of SI and S2. However, 
the emission is elongated in the direction of both mass 
structures. In the south, the ICM appears elongated in 
the direction of SI, while in the north the emission is 
elongated along the line connecting the NE and SW sub¬ 
clusters and then appears to curve in the direction of S2. 
We note that if SI and S2 have already merged with the 
NE and SW subclusters, the dark matter would have de¬ 
coupled from the gas, and therefore we do not necessarily 
expect a spatial overlap between the dark matter and gas 
components. 

Interestingly, the unsharp-masked image enchances a 
small cavity with a diameter ^ 50 kpc NW of the core of 
the NE subcluster. We examined the significance of the 
cavity from the azimuthal surface brightness profile in 
an annulus around the cluster center. The annulus was 
divided in 14 partial annuli with equal opening angles. 
The partial annuli and the azimuthal surface brightness 
profile are shown in Eigure 9. The azimuthal surface 
brightness profile is lowest in 4 partial annuli located 
NW of the NE core. The largest dip is in the partial 
annulus that crosses the middle of the X-ray cavity seen 
in the unsharp-masked image; the opening angles of this 
partial annulus are 51 — 77 degrees. No radio emission 
fills the X-ray cavity, but the cavity was likely inflated 
by the AGN hosted by the NE brightest cluster galaxy 
(see Section 8). 

8 . RADIO RESULTS 

The JVLA 1 — 2 GHz images reveal several compact 
sources in the cluster region (Eigure 10). Two of these 
sources are associated with cD galaxies in the NE and 
SW subclusters. These sources are also detected in the 
GMRT 610 MHz image (Eigure 11). These two point-like 
AGN have 1.5 GHz integrated flux densities^^ of 1.47 ± 
0.08 (NE) and 0.27 ± 0.03 (SW) mJy. At 610 MHz we 
measure flux densities of 3.24±0.34 (NE) and 0.33 ±0.07 
(SW) mJy for these sources. 

We also And diffuse extended emission in the clus¬ 
ter. The diffuse emission reveals itself in the JVLA im¬ 
age as an increase in the “noise” in the general cluster 
area. This diffuse emission is better visible in our low- 
resolution tapered image with emission from compact 
sources subtracted (Eigure 10). The diffuse emission has 
a elongated shape measuring about 120" x 45" (0.65 by 
0.24 Mpc) and is oriented along a NE-SW axis, follow¬ 
ing the overall distribution of the X-ray emission. We 
also And evidence for this diffuse emission in the GMRT 
610 MHz image, although less clearly than in the JVLA 
1.5 GHz image (Eigure 11). In the low-resolution tapered 
GMRT image there is again evidence for diffuse emission, 
but the peak flux is only at a level of 3crrms- 

Eor the integrated flux of the diffuse emission we mea¬ 
sure 1.58 ± 0.13 mJy (at 1.5 GHz) in an ellipse with 
radii of V by 0.5' orientated with a position angle of 45° 
(following the overall brightness distribution). Erom the 

Smoothing with Gaussians of larger widths does not reveal 
additional substructure. 

The quoted flux errors include a 5% uncertainty for the abso¬ 
lute flux scale bootstrapped from the primary calibrator sources. 
For the GMRT observations, we assume an uncertainty in the ab¬ 
solute flux scale of 10%. 
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Figure 5. /3-model (left) and double /3-model (right) fits to the surface brightness profiles of the NE sectors with position angles of 
45° — 115° (top), 115° — 155° (middle), and 155° — 225° (bottom). The position angles are measured from W, in a counterclockwise 
direction. For all the profiles, a double d-model constitutes a better fit than a single /3-model, with confidence levels of 4.4(7 for the top 
profile, 3.3(7 for the middle profile, and 6.4(7 for the bottom profile. The individual /3-models are shown with dashed lines. For clarity, all 
profiles shown have been regrouped in the plots to have 200 counts per bin. The best-fitting model parameters are listed on the plots; radii 
units are arcmin, and surface brightness units are photon cm“^ s“^ arcmin”^. Fixed parameters are shown with a superscripted f The 
bottom panel shows the residuals of the fit. 


GMRT tapered image we estimate an integrated flux of 
6.8 ± 3.0 mJy for the diffuse emission in the same area. 

Taking the two flux measurements at 1.5 and 0.61 GHz, 
we obtain a spectral index of 0^1500 = —1.6 ± 0.5 for the 
diffuse emission. Based on the above results, we compute 
a radio power of Pi .4 ghz = (1.06 ±0.09) x 10^^ W Hz“^, 
scaling with the spectral index of a = — 1 . 6 . 

VLITE shows emission co-incident with the NE sub¬ 
cluster where the higher resolution 1 — 2 GHz images 
(Eigure 10) reveal the compact sources and the brightest 


portion of the diffuse emission. We fit the 320-360 MHz 
VLITE emission with a single Gaussian component to 
measure the integrated flux. We measure a total flux of 
17 ± 6 mJy where we have included an 8 % uncertainty 
in the absolute flux scale. The VLITE emission contains 
both the diffuse component and the compact emission 
from the two point sources associated with the NE and 
SW clusters. We use the VLA and GMRT flux measure¬ 
ments to determine the spectral index of the two com¬ 
pact sources, assume that spectral index extends to the 
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Figure 6. Top: Sector used to model the surface brightness of the 
SW subcluster. Annuli are drawn only to guide the eye and do not 
reflect the actual bin size used for the surface brightness profiles. 
Bottom: /3-model fit to the surface brightness profile of the SW 
subcluster. For clarity, the profile shown in this plot was binned 
to have 200 counts per bin. The best-fitting model parameters 
are listed; radii units are arcmin, and surface brightness units are 
photon cm“^ s“^ arcmin”^. Fixed parameters are shown with a 
superscripted f- The bottom panel shows the residuals of the fit. 

central VLITE frequency, and estimate the contribution 
to the total VLITE emission from the compact sources. 
We subtract that from the VLITE flux and get an es¬ 
timate of the diffuse component detected by VLITE of 
11 ± 4 mJy. Comparing the VLITE flux to the JVLA 
flux, we calculate a spectral index of (afloo = —1.3 ± 0.3. 
The results are consistent with the spectral index esti¬ 
mate from the GMRT data. Deeper low-frequency data 
are required to better constrain the spectral properties 
of the diffuse emission. 

In the JVLA 1 — 2 GHz data we did not detect dif¬ 
fuse polarized emission in the cluster in Stokes Q and 
U images we made. Halos are generally unpolarized at 
the few percent level or less so this result is not surpris¬ 
ing (e.g., Eeretti et al. 2012). We note however, that a 
proper search for diffuse emission, taking full advantage 
of the large bandwidth, would require Earaday Rotation 
Measure Synthesis. 

9. GAS-DM DECOUPLING 

The combined X-ray and lensing analysis of Jauzac 
et al. (2015) determined that the gas component of the 
SW subcluster has decoupled from the DM component 



R [arcmin] 

Figure 7. Broken power-law model fit to the surface brightness 
profile of the SW subcluster in a subsector with position angles 
280° — 340°, measured from the W in counterclockwise direc¬ 
tion. For clarity, the profile shown in this plot was binned to 
have 80 counts per bin. The best-fitting model parameters are 
listed; radii units are arcmin, and surface brightness units are pho¬ 
ton cm“^ s“^ arcmin”^. Fixed parameters are shown with a su¬ 
perscripted f. 
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Figure 8. Unsharp-masked X-ray image of MACS J0416.1-2403, 
created by dividing the difference of two 0.5 — 4 keV fluxed im¬ 
ages convolved with Gaussians of widths 4" and 10'' by their sum. 
Point sources have not been subtracted from this image, and there 
are some residuals surrounding them. The positions of the two 
less massive mass structures identified by Jauzac et al. (2015) are 
marked as SI and S2 with circles of diameters 100 kpc, as also done 
by Jauzac et al. (2015). Black crosses mark the centers of the DM 
halos of the two main subclusters (M. Jauzac, priv. comm.). The 
dashed arc shows the position of the density discontinuity detected 
near the SW subcluster. The location of the X-ray cavity is also 
marked. 
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Figure 9. Left: Unsharp-masked image as in Figure 8, with overlaid partial annuli used to evaluate the azimuthal surface brighness profile 
around the NE core. Right: Azimuthal surface brightness profile around the NE core. The largest dip in the profile is in the direction of 
the cavity seen in the unsharp-masked image, between opening angles of 51 — 77 degrees. Angles are measured counterclockwise, starting 
from the W. 
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Figure 10. Left: JVLA 1 — 2 GHz high-resolution image showing the compact sources in the cluster region. Right: JVLA 1 — 2 GHz image 
of the radio halo, with compact sources subtracted. Chandra contours are overlaid on both images. The contours are based on an exposure- 
corrected, vignetting-corrected image that was smoothed with a Gaussian kernel of width 4"; they are drawn at [0.6,1.3, 2.0, 2.7,...] x 10“^ 
photons cm“^ s“^. Green crosses mark the centers of the DM halos. Dashed blue lines show the —3a radio contours (these are only visible 
in the high-resolution radio image). The beam size is shown in the bottom left corner of the images. 


and lags behind the subcluster’s galaxies as it travels 
towards the NE subcluster. In the NE subcluster, on 
the other hand, the DM and gas components were de¬ 
termined by Jauzac et al. (2015) to spatially overlap. 
Offsets between the DM and gas components of merg- 

An analysis of the DM-gas offset in MAGS J0416.1-2403 
was carried out more recently by Harvey et al. (2015). For 
the SW subcluster, their DM peak is offset by about 100 kpc 
(~ 20") from the DM peak determined by previous analy¬ 
ses (e.g., Jauzac et al. 2015; Grillo et al. 2015, models at 
https: / / archive.stsci.edu / pub/hlsp / frontier / macs0416/models/). 
We do not understand the cause of the offset, but choose to use 
the DM peaks determined by Jauzac et al. (2015), because their 
locations are consistent within 2" with the locations of the DM 
peaks determined by other authors (e.g., Grillo et al. 2015). We 


ing galaxy clusters set crucial constraints on the merger 
state. In particular, significant offsets indicate a post¬ 
merging system, while a lack of offsets supports a pre¬ 
merging scenario. Below, we discuss the DM-gas offsets 
in light of the deeper Chandra data. 

The centers of the DM halos of the two subclusters 
(M. Jauzac, priv. comm.) are marked on Eigure 12, 
which shows an HST ACS/E814W image of the cluster, 

also point out that the peak of the galaxy distribution chosen by 
Harvey et al. (2015) for the SW subcluster in MAGS J0416.1-2403 
is strongly biased by a foreground galaxy {z = O.lO^Q'Jg; Jouvel 
et al. 2014) at RA = 04 : 16 : 06.833 and DEG = -24 : 05 : 08.40 
that was not excluded from their analysis. 
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Figure 11. GMRT 610 MHz high-resolution image, showing the 
compact radio sources. Overlaid in blue are GMRT 610 MHz con¬ 
tours of the diffuse radio emission, drawn at 1.1 X 10^ /iJy/beam. 
JVLA 1 — 2 GHz radio contours showing diffuse radio emission are 
drawn in magenta at [60, 80,100,120,140,160] p,Jy/beam. 


with overlaid Chandra and JVLA contours. The centers 
of the X-ray cores were defined as the peaks in the X- 
ray emission, and were determined from the 0.5 — 4 keV 
Chandra surface brightness map of the cluster, convolved 
with a Gaussian of a = 4". We also considered the un¬ 
certainties on the surface brightness distribution, and de¬ 
fined the uncertainties on the X-ray peaks as the regions 
around the cores in which the X-ray brightness is consis¬ 
tent (within the propagated Poisson errors on the counts 
image) with the brightness of the corresponding peak. 
We note that our approach only provides a lower limit on 
the extent of the peak regions, because the background 
was included in the surface brightness map used for this 
analysis. The uncertainties on the DM centers are < 1", 
significantly smaller than the uncertainties on the cen¬ 
ters of the X-ray cores. In Figure 13, we compare the 
position of the DM centers with that of the X-ray peak 
regions. Rather than confirming the previously-reported 
offset between the gas and DM components of the SW 
sub cluster. Figure 13 shows that both X-ray peaks are, 
within the uncertainties, co-located with the DM cen¬ 
ters. We speculate that the discrepancy between the 
results presented herein and those reported by Jauzac 
et al. (2015) is caused by a higher noise level in the sig¬ 
nificantly shallower X-ray data used in previous studies; 
our data is ~ 6 times deeper, which is equivalent to a 
reduction in noise by a factor of ~ 2.5. 

10. MERGER SGENARIO 

Jauzac et al. (2015) showed that the galaxy redshift 
distribution in MACS J0416.1-2403 is bimodal relative 
to the mean redshift, with the SW subcluster (mean 
redshift 0.3966) moving towards the observer and the 
NE subcluster (mean redshift 0.3990) moving away from 
the observer. Based on these observations, Jauzac et al. 
(2015) suggested two possible merger scenarios: 


• MACS J0416.1-2403 is a pre-merging system, in 
which the SW subcluster comes from behind the 
NE sub cluster, and is now seen near first core pas¬ 
sage; 

• MACS J0416.1-2403 is a post-merging system, in 
which the SW cluster approached the NE subclus¬ 
ter from above, and is now seen near its second core 
passage. 

In this section, we try to distinguish between the two 
scenarios based on our X-ray and radio findings. 

10.1. Cooler gas in the NE? 

Approximately 300 kpc NE of the NE core, the 
temperature maps show a region (region #10; T = 
8.411q gQ keV) that is cooler than the neighboring region 

closer to the cluster center (region #8; T = 12.18l^‘9j). 
The cool gas in this region could be pre-shock gas ahead 
of a shock front. We examined the surface brightness pro¬ 
file across the cooler region, but could not confirm a sur¬ 
face brightness discontinuity with a confidence > 90%. 
A density jump could be masked by projection effects. 
Nonetheless, while the temperature difference between 
regions #8 and #10 is significant at 90% confidence (but 
not at 2<j) and there is a pseudo-pressure jump between 
the two regions, the lack of a clear density jump does not 
allow us to confirm a shock front NE of the NE core. 

10.2. No gas-DM offset 

In head-on binary mergers, the essentially collisionless 
DM travels ahead of the collisional gas. There are several 
scenarios that could explain the spatial overlap of the DM 
and gas in MACS J0416.1-2403: 

(i) MACS J0416.1-2403 is not a head-on merger that 
progresses outside the plane of the sky, and the NE 
and SW subclusters are seen before first core pas¬ 
sage and have yet to interact strongly with each 
other. 

(ii) The NE and SW subclusters are seen after first core 
passage, and the DM and gas overlap only in projec¬ 
tion. This requires either that one of the subclusters 
was essentially undisturbed by the merger and our 
line of sight aligns with the trajectory of the other 
subcluster, or that the subclusters’ trajectories are 
parallel and both are aligned with our line of sight. 

(iii) The NE and SW subclusters are seen after first core 
passage, and the gas of the SW subcluster slingshot 
and caught up with the DM halo, after they had 
separated near pericenter. 

The radial velocity difference between the BCGs of the 
NE and SW subclusters was reported by Jauzac et al. 
(2015) to be ~ 800 km s“^. The sound speed in a cluster 
with a temperature of 10 keV is ^ 1300 km s“^, and 
typical collision velocities are 1 — 2 times the sound speed. 
Therefore, the radial velocity difference between the two 
BCGs in MACS J0416.1-2403 seems rather low. The 
velocity could be explained if the subclusters are seen 
well before or after first core passage, possibly near the 
turnaround point. 
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Figure 12. HST ACS/F814W optical image of MACS J0416, with overlaid Chandra (left) and JVLA (right) contours. X-ray contours 
are the same as in Figure 10. Radio contours are based on the JVLA image with the emission from compact sources subtracted using a 
Gaussian uv-taper of 20". The JVLA contours are drawn at [60,80,100,120,140,160] /iJ/beam. The centers of the DM halos of the NE 
and SW subclusters are marked with red crosses. 
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Figure 13. Chandra 0.5 — 4 keV surface brightness maps of the NE subcluster center (left) and of the SW subcluster center (right). The 
overexposed regions (in brown) are the peak regions of the two subcluster cores. The position of the DM centers are marked with white 
crosses. The white ‘X’ indicates the location of the X-ray peak of the SW subcluster, as it was identified from shallower Chandra data by 
Jauzac et al. (2015). 


10.3. Substructure in the NE core 

The surface brightness profile of the NE subcluster can¬ 
not be described by a single /d-model, but instead by a 
double ;d-niodel composed of a very compact core and 
a more extended gas halo. Such a profile is typically 
used to model cool core clusters (e.g. Jones & Forman 
1984; Ota et al. 2013), but it can also describe the sur¬ 
face brightness of merging clusters seen in projection (e.g. 
Zuhone 2009; Machacek et al. 2010). The core of the NE 
sub cluster is very hot and has a long cooling time, which 
disfavors the possibility that the NE subcluster hosts a 
cool core. Instead, the observed substructure is more 
likely the result of a merger event. 

The most straight-forward scenario would be that the 
NE and SW sub clusters have already merged, and while 


the core of the SW subcluster was destroyed in the colli¬ 
sion, as evidenced by the subcluster’s flat surface bright¬ 
ness, the core of the NE subcluster survived. Cool cores 
that survive a recent merger event preserve part of their 
cold gas. However, we found no evidence of cool gas in 
the NE core. Therefore, the NE core is unlikley to be 
the remnant of a recent interaction between the NE and 
SW subclusters. 

An alternative scenario is that the NE subcluster is 
merging with the dark matter halo S2 8. However, S2 is 
not associated with a concentration of galaxies, and thus 
is an unlikely candidate for a group of galaxies; instead, 
S2 is probably a cosmic filament seen in projection along 
the line of sight (Jauzac et al. 2015). No other mass 
concentration that the NE subcluster might be currently 
merging with has been detected. 












16 


Ogrean et al. 


One other possibility, postulated by numerical simu¬ 
lations of Poole et al. (2008), is that the NE subcluster 
underwent a merger or a series of mergers over ^ 1 Gyr 
ago, and while its core has already relaxed back to a 
compact state, it has not yet recooled. In this scenario, 
the merger could not have been with the SW subcluster, 
because the two subclusters are involved in an ongoing 
merger. Instead, the NE subcluster needs to have inter¬ 
acted several Gyr ago with one or more clusters that are 
no longer directly detectable in the X-ray maps, in the 
lensing maps, or in the redshift distribution of the NE 
subcluster. 

Alternatively, past mergers could also have heated the 
cool core without destroying it, by inducing sloshing in 
the central galaxy. The kinetic energy of the sloshing 
central galaxy would then have been dissipated as heat 
between successive mergers. 

10.4. The AGN and the X-ray eavity in the NE eore 

The collision velocity between the NE and SW subclus¬ 
ter is unlikely to be highly supersonic. If the cluster is 
a post-merging system, a cavity could have had enough 
time to form since first core passage. However, given 
the short timescale and the intense turbulence following 
the moment of first core passage, forming a new cavity 
would require a strong AGN outburst. The expectation 
would then be to observe radio emission associated with 
the X-ray cavity. However, this radio emission is not ob¬ 
served. Instead, the presence of an X-ray cavity NW of 
the NE core supports a pre-merger scenario. The cavity 
was most likely inflated by a recent weak outburst of the 
AGN detected in the BGG of the NE subcluster. 

10.5. SW density diseontinuity 

The density discontinuity detected S of the SW sub¬ 
cluster is located along the line connecting the NE sub¬ 
cluster, the SW subcluster, and SI. The discontinuity 
also appears coincident with the position of SI, but this 
could be a projection effect. There are two scenarios that 
would explain the origin of the density discontinuity: the 
merger of the SW subcluster with SI, or the merger of 
the SW subcluster with the NE subcluster. 

If the discontinuity were triggered by a merger between 
the two main sub clusters, then the merger must have al¬ 
ready progressed past the moment of the first core pas¬ 
sage. Eurthermore, because the NE core is very com¬ 
pact and hosts an X-ray cavity, it is unlikely that it was 
significantly affected by a merger with the SW subclus¬ 
ter. Therefore, taking into account the constraints set 
by the lack of an offset between the gas and DM peaks 
ojf both subclusters, it is the SW subcluster whose tra¬ 
jectory must be aligned with our line of sight. In this 
geometry, a shock front would be expected ahead of the 
SW subcluster (with respect to its direction of motion), 
i.e. our line of sight should be perpendicular to the 3D 
shock “cap”. Instead, however, the surface brightness 
discontinuity we detect is seen in projection only S-SW 
of the SW subcluster. If the detected density discontinu¬ 
ity is the appearance of the shock in projection, then we 
would expect to observe it over a larger angular opening. 
Most likely, the density discontinuity is caused by the 
interaction of the SW subcluster with SI. 

In conclusion, while we cannot completely elim¬ 
inate the possibility of MACS J0416.1-2403 being 


a post-merging cluster, we believe the sum of our 
findings favors a pre-merging scenario. 

11. A NEWLY-DISCOVERED RADIO HALO 

We classify the diffuse radio emission in the cluster 
as a halo, since the emission has low-surface brightness, 
large extent, and is centrally located. The halo has a 
somewhat smaller size than other radio halos in merging 
clusters, which have typical sizes of 1 — 1.5 Mpc. (e.g., 
Eeretti et al. 2012). The size is more similar to that of 
radio mini-halos found in cool-core clusters (e.g., Giacin- 
tucci et al. 2014). However, an important difference from 
mini-halos is that neither the NE nor the SW subclusters 
contains a cool core. In addition, as shown in Eigure 12, 
the halo seems to be associated with both subclusters. 

The radio emission around the northern subcluster 
contains about a factor of 2 more flux than the emis¬ 
sion around the southern one. An interesting question 
is whether we observe a single radio halo, caused by the 
merger between the NE and the SW sub clusters, or two 
individual halos belonging to the two subclusters. In the 
latter case, these halos must have originated from pre¬ 
vious merger events within two sub clusters themselves. 
This would make MAGS J0416.1-2403 similar to the case 
of the radio halo pair in Abell 399 and Abell 401 (Mur- 
gia et al. 2010), with the difference that Abell 399/401 
seems to be at a much earlier pre-merger stage. 

Eor radio halos, a correlation is observed between the 
cluster’s X-ray luminosity and the radio halo power (the 
I/x — P correlation), with the most X-ray luminous clus¬ 
ters corresponding to the most powerful radio halos (e.g., 
Liang et al. 2000b; EnBlin & Rottgering 2002; Gassano 
et al. 2006). The X-ray luminosity is used as a proxy of 
cluster mass. The distribution in the Lx — P plane is bi- 
modal when clusters without halos are included. Merging 
clusters with radio halos follow the Lx — P correlation, 
while the upper limits on the radio power of other clus¬ 
ters are well below the correlation (Brunetti et al. 2007). 

A correlation is also observed between the cluster’s 
integrated Sunyaev-Zel’dovich (SZ) effect (i.e., the in¬ 
tegrated Compton Tsz parameter) and the radio halo 
power (e.g., Basu 2012, 2013). Tsz is proportional to the 
volume integral of the pressure. The SZ signal should be 
less affected by the cluster dynamical state and therefore 
is expected to be a better indicator of the cluster’s mass. 
Basu (2012, 2013) and Gassano et al. (2013) showed that 
there is good agreement between the Tsz — P scaling re¬ 
lation and the scaling relations based on X-ray data. 

In Eigure 14, we place MACS J0416.1-2403 on the 
Lx — P and Tsz — P diagrams, with the values for 
the other clusters taken from Gassano et al. (2013). 
The Tsz, 500 value for the cluster was recently calcu¬ 
lated by Planck Collaboration et al. (2015), which found 
^sz, 500 = (0.71 ± 0.24) X 10“^ Mpc^. Specifically, this 
is the value of Tsz, 500 obtained from the union cata¬ 
log based on the MMEl^^ detection algorithm. We also 
computed the value of Tsz, 500 using the gNEW model 
fit to Bolocam presented by Czakon et al. (2014) and ob¬ 
tained Tsz, 500 = (1.49 ± 0.28) x 10“^ Mpc^ The differ¬ 
ence between the Planck and Bolocam measurements is 
not understood, but is likely due to some combination of 
random noise fluctuations, deviations from the assumed 

Matched multi-filter algorithm, implementation 1. 
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Figure 14. Lx — P and Isz ~ P diagrams. Clusters with only upper limits on the radio halo power are shown in teal-colored circles, 
while clusters with detected halos are shown in dark red squares. On both diagrams, MACS J0416.1-2403 is represented with a star. In 
the Isz — P diagram, the orange star corresponds to the Isz, 500 value reported by Planck Collaboration et al. (2015), while the blue star 
corresponds to the Isz, 500 value calculated from the gNFW fit to Bolocam presented in Czakon et al. (2014). Error bars are drawn for 
all points, but some are very small and barely visible. Values for all clusters with the exception of MACS J0416.1-2403 are taken from 
Cassano et al. (2013). An interactive version of the diagrams is available at http://hea-www.cfa.harvard.edu/~gogrean/Lx-P.html. 


pressure profile used to extract Isz, 500 (particularly at 
large radii), coutauiiuatiou froui dust or other astrououi- 
ical siguals, aud calibratiou uucertaiuties. Furtheruiore, 
we uote that this cluster is uot detected by PowellSuakes 
- oue of the three Plauck algorithuis - aud the values of 
^sz, 500 obtaiued froui the other two algorithuis (MMFl 
aud MMF3^^) are siguificautly differeut, perhaps iudi- 
catiug that uieasureuieuts of Isz, 500 towards this cluster 
coutaiu larger thau typical systeuiatic uucertaiuties. lu 
Figure 14, we plot both the Plauck aud Bolocaui values 

of ^SZ, 500- 

luterestiugly, the radio power falls ou the low eud of 
the Lx — P correlatiou, which uieaus that the radio 
halo iu MACS J0416.1-2403 is souiewhat less luuiiuous 
thau expected froui the X-ray luuiiuosity of the cluster. 
However, usiug the Plauck uieasureuieut of Fsz, soo, we 
hud that the halo power is cousisteut with the expected 
^sz — P correlatiou, iudicatiug the halo is uot uuderlu- 
uiiuous for the cluster uiass. While the Bolocaui value 
of hsz, 500 does iuiply lower thau expected halo power, 
the uouiiual Isz — P relatiou used for couiparisou is de¬ 
rived solely froui Plauck SZ uieasureuieuts. We therefore 
cousider our iuterpretatiou based ou the Plauck value 
of Tsz, 500 to be uiore robust. Giveu this iuterpreta¬ 
tiou, MACS J0416.1-2403 could be just au outlier iu the 
Lx — P diagraui, based ou the fact that there is siguifi- 
caut scatter iu the correlatious. 

It is also possible that the X-ray luuiiuosity of 
MACS J0416.1-2403 is uot a good iudicator of the 
cluster’s uiass, aud the luuiiuosity is boosted by the 
uierger. A boost iu X-ray luuiiuosity is predicted 
for approxiuiately oue souud-crossiug tiuie, which for 
MACS J0416.1-2403 is ^ 1 Gyr (Ricker & Saraziu 
2001). This boost would affect other clusters ou the 
correlatiou as well, but it could affect MACS J0416.1- 
2403 uiore siuce we are seeiug a cluster close to first 
core passage. The uiass of the cluster iu L 500 is 
(7.0 ± 1.3) X 10^^ Mq (Uuietsu et al. 2014). Based ou 


the scaliug relatious of Pratt et al. (2009), the uiass 
correspouds to a 0.1 — 2.4 keV cluster luuiiuosity of 
(1.1 ± 0.4) X 10^^ erg s“^, which is iu agreeuieut with 
the luuiiuosity calculated froui the Chandra data. There¬ 
fore, based ou the M 500 — L 500 scaliug relatiou, we have 
uo iudicatiou that the cluster luuiiuosity is higher thau 
expected for a cluster of this uiass. 

Auother possibility is that we are seeiug au ultra-steep 
spectruui radio halo {a < —1.5, Bruuetti et al. 2008). Iu 
this case, the radio halo is ouly uuderluuiiuous at higher 
frequeucies (i.e., above a GHz), while the iutegrated flux 
deusity rapidly iucreases at lower frequeucies. The 320 — 
360 MHz VLITE aud 610 MHz GMRT data preseuted 
iu Sectiou 8 are uot sufficieut to coufirui, or rule out, 
the preseuce of au ultra-steep spectruui radio halo. Low- 
frequeucy GMRT {u < 610 MHz) aud/or deep P-baud 
JVLA observatious are required to settle this case. 

A third possibility is that we could be wituessiug the 
first stages of the foruiatiou of this radio halo. We do uot 
thiuk that the radio halo is fadiug, because the NE aud 
SW subclusters are still iu the process of uiergiug aud 
uiost likely they have uot yet uudergoue core passage 
(see Sectiou 10). Therefore, wheu the uierger progresses 
further, the radio halo could brighteu, iucrease iu size, 
aud uiove outo the correlatiou. 

12. SUMMARY 

The HST Eroutier Eield cluster MACS J0416.1-2403 
{z = 0.396) is a uiergiug systeui that cousists of two uiaiu 
subclusters aud two additioual suialler uiass structures 
(Jauzac et al. 2015). Here, we preseuted results froui 
deep Chandra JVLA, aud GMRT observatious of the 
cluster. The uiaiu aiuis of our aualysis were to ideutify 
siguatures of uierger activity iu the ICM, coustraiu the 
uierger sceuario, aud detect aud characterize diffuse radio 
sources. Below is a suuiuiary of our uiaiu fiudiugs: 

• The cluster has au average teuiperature of 
10.06lo!49 keV, au average uietallicity of 

0.241qq4 Zq, aud a 0.1 — 2.4 keV luuiiuosity 
of Lx = (9.14 ± 0.09) X 10^^ erg s-k 


Matched multi-filter algorithm, implementation 3. 
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• The NE subcluster has a very compact core and 
a nearby X-ray cavity, but its temperature is very 
high (T ^ 10 keV), we find no evidence of multi¬ 
phase gas, and its cooling time is > 3 Gyr. The 
presence of a compact core and of an X-ray cavity 
indicates that the NE subcluster was not disrupted 
by a recent major merger event. 

• S-SW of the SW subcluster, there is a weak den¬ 

sity discontinuity that is best-fitted by a density 
compression of l.SGloJg* discontinuity is lo¬ 

cated along the line connecting the NE subcluster, 
the SW subcluster, and the mass structure SI re¬ 
ported by Jauzac et al. (2015). Most likely, the 
discontinuity was caused by a collision between the 
SW subcluster and SI, rather than by a collision 
between the NE and SW sub clusters. 

• The DM and gas components of the NE and SW 
subclusters are well-aligned, which favors a sce¬ 
nario in which the two subclusters are seen before 
first core passage. 

• MACS J0416.1-2403 has a radio halo with a 
1.4 GHz power of (1.06 ± 0.09) x 10^"^ W Hz“^, 
which is somewhat less luminous than predicted by 
the Lx — P correlation. However, the halo aligns 
well on the Fsz, 500 ~ P correlation, indicating that 
it is not underluminous for the cluster mass. We 
could be observing the cluster at a point in the 
merger event when the X-ray luminosity is signifi¬ 
cantly boosted. Alternatively, we could be observ¬ 
ing the birth of a giant halo, or a very rare ultra¬ 
steep halo. 
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